MAXIMIZING A COMBINATORIAL EXPRESSION ARISING 
FROM CROWD ESTIMATION 



DONALD M. DAVIS 



Abstract. We determine, within 1, the value of N for which 
S ( S i) (jv) (n-i) ( ■ ) achieves its maximum value. Here s% and S2 

i 

are fixed integers. This problem arises in studying the most likely 
value of |j4U-BUC| if A and C are disjoint sets of cardinality s\, and 
\B\ = S2- Attempting to remove the 1 unit of indeterminacy leads 
to interesting conjectures about a family of rational functions. 

1. Introduction 

The question considered here arises from problems involving estimating sizes of 
crowds. You count sizes of certain subsets and want to estimate the size of the union. 

The case considered here involves three sets A, B, and C with the property that 
A n C = 0, but B may intersect the other sets. Suppose also that \A\ = s\, \B\ = s 2 , 
and \C\ = S3. What is the most likely value for \A U B U C|? This question was 
suggested to the author by Fred Cohen, along with the mathematical model which 
we now present. 

The assumption being made is that any choice of i people in A D B and j other 
people in B n C is equally likely, for all % and j. For example, it is equally likely that 
\A n B\ = 1 with that person a specified person of A and a specified person of B or 
that \A fl B\ = 3 with those people a specified subset of A and a specified subset of 
B. This can be formulated in the following way. 

Suppose A, B, and C are disjoint sets with \A\ = si, \B\ = s 2 , and |C| = s 3 . The 
sample space consists of all 4-tuples 

{A^B^B^O,) c (A 
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such that I Ai| = \Bi\, \B 2 \ = |Caj , and B\ D B 2 = 0. Thus A\ and -Bi correspond to 
Ad B in the earlier formulation, and we have 

\A U B U C\ = si + s 2 + s 3 - |Bi| - |S 2 |- 

We assume that each element of the sample space is equally likely. Let Eij denote 
the event that | Bi \ — i and \B 2 \ — j. Then 

i-^ (?)(?) (7*) (?)■ 

If E N is the event that \B X \ + \B 2 \ = N, i.e., that \Al>Bl>C\ = s x + s 2 + s 3 - N, then 

i+j=N 

Hence the most likely value of \A U B U C\ is s\ + s 2 + S3 — N, where N maximizes 

E(")(?)(^C^ = E(?)S)(A)(7)- 

i i 

We focus attention primarily on the case in which s± — S3. In this case we obtain 
in Corollary 11.31 a simple formula for the maximizing iV within 1 , and in 11.101 a much- 
less-tractable formula which removes the indeterminacy. In Section [31 we attempt to 
obtain a more useful approximation to I1.10[ and in doing so we notice fascinating 
patterns in a family of rational functions, but can only conjecture that these patterns 
persist. See Table |3] and Conjecture 13.91 In Section @J we consider the general case 
when Si and S3 need not be equal. Our results there are somewhat similar, but not 
so complete. 

Our main theorem is 

Theorem 1.1. Let f Sl , S2 (N) := (£) E (?) G£j (?) f or integer values of N. For 

i 

each Si and s 2 , there is an integer, which we denote by g{si,s 2 ), such that / sl) s 2 (iV) 
is an increasing function of N for N < g(si, s 2 ), and a decreasing function of N for 
N > g{ s ii s 2)- Moreover, 



(1.2) 9(8!, 8 2 ) 

with 5 = or 1. 



2si + s 2 + § - \/4:sj + 4si + (s 2 + \f 



+ 5 
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Corollary 1.3. The maximum value of ^2 (7) (jy) Gv-0 { N ) 



occurs when 




with 5 = or 1. The most likely value of \AU BUC\ in the situation discussed above, 



It is conceivable that / S1)S2 might achieve equal maxima at both N and N + 1. In 
such a case, we accept either as an allowable value of g(s\, s 2 ). 

To illustrate the efficacy of our formula, we consider the typical case s\ = 15, and 
tabulate in Table [TJ the actual values of <?(15, s 2 ) for all S2, and in Table [2] the five 
values of S2 for which 5 = 1 in (ll.2p . Note how in these five cases the expression 
whose integer part appears in (11. 2p falls slightly short of the required value. 

The following proposition generalizes the beginning and end of Table [TJ Theorem 
11.11 is true with 5 = in these cases. Note also that the case d = of part (b) of 
Proposition 11.41 shows that if B is much larger than A and C, then the most likely 
occurrence is that both A and C are contained in B. 

Proposition 1.4. 

a. If s 2 < |(V 8s i + 9 - 1); then g(si, s 2 ) = s 2 . 

b. For a// Si, s 2 , we have g(si, s 2 ) < 2si. For < d < 4, 




3 
2 



5. 



(1.5) 



»(si, s 2 ) >2 Sl -d iff s 2 > -^{s\ + si) 



d+2 
2 
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Table 1. Values of g(15, s 2 ) 



$2 


#(15, s 2 J 


f1 fil 


S 2 




o 1 

s 2 — 1 


[11, 12J 


c 9 
S 2 — Z 


[id, loj 


s 2 — 


r 1 1 7i 

[ID, 1/j 


-1 
S 2 — 4 


1 & 
lo 


1 Q 


ri q oil 
[iy, zij 


1 A 


[00 OQl 
[ZZ, ZoJ 


1 ^ 
10 


fo/1 9fil 
[Z4, ZOJ 


1 

10 


[97 9Ql 

[z I , zyj 


1 7 
1 ( 


TQn qqI 
[oil, ooj 


lo 


[o4, O (J 


1 Q 

iy 


Too /|Ol 


9n 
zu 


Li Q /I Si 
[4d, 4SJ 


91 
Zl 


[4y, OOJ 


99 

zz 


[00, 04tJ 


9^ 

zo 


65 76 


24 


[77, 92] 


25 


[93,117] 


26 


[118, 157] 


27 


[158,238] 


28 


[239, 478] 


29 


> 479 


30 



Table 2. Cases in which equality does not hold in (11.21) 
when Si = 15 and 5 = 



S2 


0(15, s 2 ) 


2si + s 2 + | - ^4s? + 4si + (s 2 + \f 


6 


6 


5.84 


10 


9 


8.78 


15 


12 


11.85 


17 


13 


12.91 


19 


14 


13.89 



Next we introduce the polynomials involved in the proof. We will usually replace 
s\ by x, both because it will occur as a variable in polynomials, and so that we can 
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use the notation Xj = x(x — 1) • • • (x — % + 1). For a nonnegative integer d, define a 
polynomial P<f(x) of degree 2d by 

(Xi) 2 (x d _i) 2 



i=0 v 1 

We will prove the following key result in Section [2j 

Lemma 1.6. When Pd+i(x) is divided by Pd(x), the quotient is j^x 2 — -^jx + |. 
Let Rd{x) denote the remainder. If f is as in \l.l\ then, if x and d are integers, 

(1.7) f XtS2 (2x-d) > f x>S2 {2x-d-l) iffs 2 > ^( x i +x )-^ + MEL. 

In Section [2], we will also prove the following result, the proof of which is less 
straightforward than that of Lemma 11.61 

Lemma 1.8. For d>5 and x > d/2, -0.5 < Rd{x) / 'Pd(x) < 0. 
Now we can prove the main theorem. 

Proof of Theorem \l.l[ By Lemma fX .81 increasing d by 1 changes Rd{x)/Pd{x) by at 
most 1/2, and clearly it decreases -^(x 2 + x) — by more than 1/2. Thus, for 
x > d/2, the RHS of (11. 7p is a decreasing function of d, at least for integer values of 
d. This, with (11.71) . implies the unimodality part of the theorem, with maximum of 
f Sl ,s 2 (N) occurring for N = 2s\ — d for the smallest integer d such that the RHS of 
(11.71) is satisfied. 
We will show that Lemma 11.81 also implies that 

(1.9) [2s x - d x \ < g( Sl , s 2 ) < [2 Sl - d 1 ] + 1, 

where di satisfies 

q — 2 (~2 i s \ dx+2 



This value is d x = -s 2 - § + J 4s 2 + As x + (s 2 + |) 2 , yielding ffL2l) . 

To prove (11.91) . write d\ = d 2 — t, with < t < 1, and d 2 an integer. Thus 
[2si — d\] = 2si — d 2 . The RHS of ( 11. 7ft is satisfied using d 2 since, using 1X781 at the 
last step, 

^2 - dT+rUi+SiJ — ^ ^TTl s i+ s iJ 2" ^ ^TTl s i+ s iJ —+pT(^- 
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Therefore g{s\, s 2 ) > 2s\ — d 2 . We will now show that the RHS of fll.7jl is not satisfied 
using d = d 2 — 2, which implies g(s\, s 2 ) < 2si — d 2 + 1, hence completing the proof. 

To see this, let h(d) = s 2 - (£r(4 + s i) ~ ¥)• Then K d i) = and di - (d 2 -2) > 
1, hence h(d 2 — 2) < — ~. Therefore when d = d 2 — 2, using Lemma 11.81 again, 
s 2 - (^(s? + si) - ^ + gg}) < 0, as desired. □ 

In terms of Rd/Pd, we give in 11.101 a precise result about whether 5 = or 1 in 
Theorem 11.11 The usefulness of this is limited by the complicated nature of Rd/Pd- 
In Section [3], we discuss a very strong conjecture regarding Rd/Pd, which, if proved, 
would make Theorem 11.101 more useful. See Theorem 13.21 The evidence for this 
conjecture leads to remarkable conjectural patterns among some rational functions. 
See Table [3] and Conjecture 13.91 

Theorem 1.10. Let 

d = 



4s? + 4si + (s 2 + |) 2 -s 2 -| 
Then $1.2\) is true with 5 = 1 iff 

(1.11) 0<^ + *i)-^-s 2 <-^S- 



Proof. Let d 1 = ^As\ + As x + (s 2 + \f - s 2 - §. Then s 2 = ^(sj + s x ) - 

If d\ is an integer, then do = d\ and the (> 0)-condition in (11. lip is not satisfied, 
and the RHS of (11.71) is not satisfied when d = d\ — 1 by an argument similar to the 
proof of 11.11 Hence g(s\, s 2 ) = 2s\ — d\, verifying Theorem 11.101 in this case. 

If d\ is not an integer, then ( 11.21) is true with 5 = 1 iff the RHS of (11. 7ft is satisfied 
using do (since 2s\ — d = [2si — d±] + 1), but the RHS of ( 11. 7ft is exactly the <-part 
of (11.111) . Note that the (> 0)-part of (11. lip is certainly satisfied in this case, since 
the middle expression in ( 11. lip equals using d\, and is a strictly decreasing function 
of d. □ 

2. Combinatorial proofs 
In this section we prove Lemma [1.61 Proposition 11.41 and Lemma [1.81 
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Proof of Lemma 11.61 Cancelling common factors in the binomial coefficients involving 
S2, we find that the LHS of (11. 7j) is equivalent to 

{x\) 2 {2x - d)\ 



(s 2 -2x + d + l)^- 



(i!) 2 ((2x -d- i)\) 2 (x - i)\(d + i-x)\ 
(x\) 2 {2x-d-l)\ 



(i\) 2 ((2x -d-1- i)\) 2 (x - i)\{d + l + i-x)V 
Cancelling (2x — d)\ and letting j = x — i, we obtain the equivalent condition 



(s 2 - 2x + d + f ) 



(x\) 2 



( x _ J )!2 (x+J _ rf)! 2 j!(rf _ j)! 



> E ' r ' )2 



(x-j)\ 2 (x+j-d-iy. 2 j\(d+i- j)v 

Multiplying both sides by (a;!) 2 , the inequality becomes 

r« O r I rll llV ( X j) 2 ( X d-j) 2 > \- (gj) 2 (Xd + l-j) 2 

and this readily yields 

(2.1) f XyS2 (2x - d) > /,, S2 (2x - d - 1) iff s 2 > + 2x - d - 1. 

Our next claim is that the leading terms of Pd{x) are given by 
(2.2) 

ft « - " (Si*"" + 3|&( M2 - M + ^ + ta »' 

The proof of (12.21) makes frequent use of 

which is true since each side is the coefficient of x d in e x ■ e x = e 2x . Then (12.31) is 
exactly the coefficient of x 2d in Pd(x). The coefficient of x 2d ~ x is 

Ei(i — 1) + (d — i)(d — i — 1) _ ^ x _ 2 d-i 
iifd-iV ~ ~ 2^ (i-2)!(d-o! - "(^yi- 

i=0 ' v ; ' i=0 

i-1 

We have used symmetry in the first step. Note that the i{i — 1) comes as ^2 %3- 

3=0 
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The next coefficient of Pd{x) is obtained similarly, but involves much more work. 
The coefficient of x 2d ~ 2 in (xi) 2 (xd^i) 2 is 

i— 1 d— i— 1 i— 1 d— i— 1 

j=0 j=0 j=0 j=0 l<ji<h<i l<h<h<d-i 

Noting that 

(i— 1 x 2 i— 1 
j-^ji-^j^^' j=i i=i 

we obtain 

" flf ~ + (« - W - ^ - 1)(^ - + 2(^) 2 + 2(^1^) 2 

j=0 j=0 

= _ i^-m^-^-i) + {i _ md _ i _ l){d _ i) 

, i 2 (i-l) 2 , (d-i) 2 (d-i-l) 2 
2 2 

Thus, using symmetry, the coefficient of x 2d ~ 2 in Y) — l — ^— is 

i=o i!(a-i)! 

E t(i-l) _,_ V- 1 1 Y" 2i— 1 

(j-2)!(d-i)! ^ (i-2)!(d-j-2)! 3 ^ (i-2)!(d-j)! " 

Next note that 

E »(i-1) \^ (j+l)(j+2) 

(i-2)!(d-i)! j!(d-2-j)! 

2 d - 2 2 d - 3 (j-l) + l 

— 9 I- 3 u \ KJ 



( d _2)! (d-3)! ^(j-l)!(d-2-j)! 

2^—1 2^ ^ 2^~^ 2^ ^ 

+ 3— — — + — — + 



(d-2)! (d-3)! (d-4)! (d-3)! 
_^il_ (8 + 8(d-2) + (d-2)(d-3)) 

od-4 

(d 2 + 3d-2). 



(d-2)! 
Also, 

E (i-2)!(d-i)! = 2 ^ (i-3)i(d-i)! + ^ (j-2)!(d-i)! = (1=2)! -2 + 3). 
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Hence the desired coefficient equals 



2<2— 4 2 d ~ 4 2 rf_2 

(d 2 + 3d-2) + — — — - ^r (d + l) 



3(d-2)! 
as asserted in (I2.2p . 
Next we claim that 



(d-2)\ K ' (d-4)! 3(d-2 

2^—3 

(3d 2 -5d + 4), 



(n A\ Pd+l{x) __ 2 2 2d_„ _i_ c[ _i_ gd(g) 

U ' " p d ( x ) d+l X 2 ^" P d (x)- 



with deg(i?d(x)) < 2d, the claim of the first part of Lemma [L6l This can be discovered 
by division of polynomials, using (12. 2p . but it is simpler just to verify that 

- ft 12 " + w^ (M2 - m + 4 >* M i - + ^ 

Combining (O) and (E3D, we obtain (O). □ 



Proof of Proposition L4_, (a). One can easily show that the bracketed expression in 
(11.21) is always less than s 2 + 1, and equals s 2 if and only if s 2 satisfies the hypothesis 
of 11.4( a). Noting that g(si,s 2 ) is > the bracketed expression of (11. 2p by Theorem 
11.11 and is < s 2 since f Sl , S2 ( s 2 + 1) = 0, part (a) follows, 
(b). The first statement is true since / sliS2 (2si + 1) = 0. 

After cancelling common factors in the numerator and denominator, one can com- 
pute that 

'0 d = 0,l 

(2x - 1) /(3(2x 2 - 2x + 1)) d = 2 

3(x-l)/(4(x 2 -x + l)) d = 3 

6{2x 3 - 9x 2 + 3x- 9)/(5(2x 4 - 8x 3 + 14x 2 - 5x + 3)) d = 4 

Thus (II. 5p follows immediately from (11.71) if d — or 1. 

Note that (ll.7p is only meaningful if 2x > d. If d = 2, then x > 2. For such x, 
—0.2 < R 2 (x) I P 2 {x) < 0. Since, for integer x, |(x 2 + x) — 2 is either an integer or an 
integer plus 1/3, an integer s 2 satisfies 

s 2 > l(x 2 + x) - 2 + iff S2 > |( x 2 + x) _ 2) 



-Rd(x) 
PAx) 
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and so (11.51) follows from ( 11. 7(1 . 

A similar argument works for d = 3 and 4. For <i = 3, we have that -£^(x 2 +x) — ^y 2 - 
is an integer plus 1/2, and for x > 2, —0.25 < R^{x)/P^{x) < 0. If d = 4 and 
x > 3, then ^-(x 2 + x) — ^j 2 - is an integer or an integer plus t with t > 0.2, while 
-0.12 < R 4 (x)/P 4 {x) < 0. If d = 4 and x = 3, ([I7D says s 2 > 2.8 - 0.55, while the 
hypothesis says s 2 > 2.8. These are, of course, equivalent. □ 

Proof of Lemma \1.8l We begin by removing common factors in Pd(x) and Pd+i(x). 
Since parity of d plays a role, we let d = 2b + e with e G {0, 1}. When considering 
R 2b+t (x), we let, for 5 G {0, 1}, 

Pzb+e+SyX) 



P 2b+e +s(x) :-- 



n 



6-l+e 
i=0 



Note that P 2 b +e+ i(x) / P 2 b+ e (x) has the same quotient as P 2 b+ e +i(x) / P 2 b+ e (x) , while 



ft 



26H 



its remainder R 2 b +e (x) satisfies R 2 b+ e (x) 

R2b+e(x) / P 2b +e(x) = R 2 b+e {x) / P 2 b+e {x) ■ 

To prove the lemma, we will prove 

(1) P 2b+t (x + b) > for x > 0, 

(2) R 2b +e(x + b) < for x > 0, and 

(3) R 2b+ e{x + b) + \P 2b+t {x + b) > for x > 0. 

We have, with q & = 1 unless z = 6, while q, & = 



^)/rii=o +e ( :r ~~ an d hence 



' 6 b-i 



i=0 j=l 

£(T)nV 

i=0 j=0 



j=6— i+l 

j? n (x+j) 2 

j=&— i+1 



6+1 b-i 

E^(T)n(^ 

i=o i=o 



j -i) 2 n (x+j) 

j=b—i+l 



e + 5 = 
e + 5=l 
e + 5 = 2. 



Part (1) is true since P 2 b+ e {x + b) is a sum of nonnegative terms including the term 

b 

Y[ (x + j) 2 , which is positive for x > 0. 

j=l-e 
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Next we consider (2) with e = 1. We compute q2b+i(x + b) = -^{x 2 — x) + \. We 
have 

^±^R 2b+1 (x + b) = ^(P 2b+2 (x + b )-q 2b+1 (x + b )P 2b+1 (x + b )) 

= Htiix-j) 2 n (*+j) 2 V*, 

i=0 \'=1 j=b-i+l ' 

where 

Fi = m(-H.-if-(2H2)( 2 r)(i(, 2 -x)+i) + iav 
= ( R-i 1 ) - ct) + *,* m ) - 2 - 2 ( (Ti) +(b-i) m ] 

+ ( 2 ^)(6+l-.) 2 -(6 + l)( 2 ^) 

( x - (b + 1 - i)f - ( 26+1 ) ( x + b - if + ( 2fe+1 ) (2(6 - i) 2 + 6 - 2i) Kb 

(b(-i) 2 -o 4 1 

Here 5^ is the Kronecker delta. 

The first term of (the last form of) and second term of when multiplied 
by the appropriate double products, cancel. Thus we obtain 

(2.5) ^flR 2b+1 ( x + b ) 



x 



2 

1 ,b- 



= e(i>-^) 2 n (-+j) 2 )m(2(6-,) 2 +&-2,) 

i=0 S'=l j=b-i+l ' 

-(ri(-+i) 2 )KT)- 

\ 7 =i / 



3 

One can easily prove that 

^( 2 Y 1 )(2(6-,) 2 + 6-2,)-K 2 ^)=0. 

i=0 

Also 2{b-i) 2 + b-2i > iff i < b+\-\^/W+l. Note also that the double products 
are increasing with i for x > 0. Thus our expression for ^ 2b ^' R 2 b+i(x _|_ ^ j s Q f the 

b 

form aj/3j with < ai < • • • < ^2 fa = 0, and fa > iff i < i . Such a sum is 

i=0 

negative. 
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The proof for (2) when e = is extremely similar. We have q2b(x+b) = jbTi^^WhT' 
Then 

i=0 j=l j=b—i+l 

where now 

Fi = ( 2b t 1 )(x-b + i r-c l ^)(2x 2 + b) 

(£) (* - (6 - 0) 2 " (?) (x + & - i? + (?) (2(6 - if - b) i < b 

( b 2 \y-Kt!) i=b. 

The rest of the argument follows exactly the same steps as in the last paragraph of 
the above proof of the case e = 1, using Y^iZo (?) (2(& — i) 2 — b) — b^Zi) = 0. 

The proof of (3) is essentially the same, except that we are subtracting 1/2 from 
q2b+e(x + b). The effect, when e = 1, is to add [b + 1)( 2 t ) to F{. The replacement 
for (J23D is 

^(R 2b+1 (x+b) + ±P 2b+1 (x + b)) 
= II (x + j) 2 )(T)(2(^-0 2 + 26-2, + l) 

i=0 \j=l j=b-i+l ' 

+(n(-+i) 2 )(T) 

\j=i / 

which is clearly positive for x > 0. The proof when e = is similar. □ 

3. Conjectures about R d (x)/P d (x) 

In Theorem 11.101 we determined the precise value of our focal function g(s±, S2) in 
terms of Rd{si)/Pd{si). In order to make this result useful, we need better informa- 
tion about the family of functions Rd(x) / Pd(x) . In this section, we present several 
conjectures about this family of functions, one of which is supported by remarkable 
patterns. See Table [3] and Conjecture 13 .91 We also discuss their implications. 

We now state the simplest of these conjectures. 



R d (x) ^ ,995(rf-2) 



Conjecture 3.1. If d > 5 and x > \{d + y/d + 2), then — > ,,,. , ,,_> . 
The implication of this conjecture is given by the following theorem. 
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Theorem 3.2. Assume Conjecture \3. II and S2 > h(y/8si + 9 — 1). Let 



d 



4sf + 4si + (s 2 + \) 2 - s 2 - | 



> 5. 



Then $1.2\) is true with 5 = 1 if 

For Sx < 38, the only cases in which (11. 2p is true with 5 = 1 which are missed by 
this theorem are (sx,s 2 ) = (6,4), (18,56), (36,16), and (38,155). The significance of 
the 0.995 in 13. II is that it is, to three decimal places, the largest number for which the 
inequality appears to be true. 



Proof of Theorem\3lM Let d = y As\ + 4sj + (s 2 + \) 2 - s 2 — §. The theorem is 
vacuously true if c?o — d because -^(sf + si) — — s 2 = 0. So we assume d is not 
an integer. Then 2s i — d = [2si — d] + 1, and so, using (11.71) . the assertion that (11.21) 
is true with 5 = 1 can be stated as 

This will follow from H 3 . 3 f) and our assumption of Conjecture 13. II once we know that 
si > ^ + ^^+2). 



Since y As\ + 4si + (s 2 + |) 2 — s 2 — § is a decreasing function of s 2l and c?o < d, it 
suffices to prove S\ > h(d + \/d~+2) if 



d = JAs\ + 4si + (| V8si + 9) 2 - i V8si + 9-1 



2sx + i - \/2si + | 



2 y ""A 1 4" 

Solving the latter equation for si yields exactly si = ~(d + y/d + 2). □ 

Extensive Maple calculation led the author to expect that, for d > 5 and x > d/2, 

R d (x) _ -(d-2) 



(3.4) 



PAx) 2x + d-2' 
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To understand how good is the approximation (13.41) . we consider the ratio of the 
two sides, using the reduced versions R and P. As this will be close to 1, we study 

R d (x)/P d (x) 



Qd(x) 



(3.5) 



-(d-2)/(2ac + d-2) 
(d - 2)P d {x) + (P d+1 (x) - q d (x)P d (x))(2x + d 
(d-2)P d (x) 



where q d (x) = -^x 1 — j^x + 1 is the quotient in (12 ,4p . We would like to prove that 
Qd{x) ~ in some sense, when x > d/2. 

Similarly to the methods in deriving (12.21) . we can show that 

(3.6) Jim o Q d (x) = {d+1) ? d _ 2y 

To see this, note that the numerator and denominator of (I3.5P are both polynomials 
of degree 2[|]. The desired limit in (13.61) is the ratio of their leading coefficients. We 
omit the details in this computation. 

We begin by considering Q25{x). It is a ratio of two polynomials of degree 24. 
Maple computes that the derivative of Q25{x) is only at x ~ 4.0409 and 60.50336. 
Moreover, Maple plots the graph of Q25(x), which turns out to look something like 
the rough sketch in the top half of Figure 13.71 This sketch is not at all to scale. We 
are particularly interested in the values for x > 13. As x increases from 13 to 60.5, 
Q 25 {x) decreases from 0.01672 to -0.006867. 




The amazing observation is that, for all odd d, Qd{x) apparently has a form very 
similar to that in the top half of Figure 13.71 with only one local maximum and one 
local minimum, which are the absolute maximum and minimum. Note that Qd{x) is 
a ratio of two polynomials of degree d — 1, yet it apparently has this simple form for 
all odd d. 
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Table 3. Max, min, and lim of Qd{x) when d is odd 



d 


X max 


Qd (^max ) 




Qd (^min) 


Qlim 


5 


0.6874697648 


3.069437405 


10.48133802 


-0.1360199330 


-0.1111111111 


7 


1.044428415 


4.216907295 


15.48679761 


-0.06528852351 


-0.05000000000 


9 


1.371314912 


5.374255987 


20.49141428 


-0.03966514187 


-0.02857142857 


11 


1.706180162 


6.540442715 


25.49465519 


-0.02724322575 


-0.01851851852 


13 


2.040065994 


7.712042960 


30.49698936 


-0.02018382313 


-0.01298701299 


15 


2.373641870 


8.886960397 


35.49873537 


-0.01574295269 


-0.009615384615 


17 


2.707177063 


10.06413875 


40.50008573 


-0.01274397486 


-0.007407407407 


19 


3.040658416 


11.24291343 


45.50115931 


-0.01060966757 


-0.005882352941 


21 


3.374105921 


12.42285407 


50.50203244 


-0.009028116216 


-0.004784688995 


23 


3.707530371 


13.60364584 


55.50275603 


-0.007817967242 


-0.003968253968 


25 


4.040938032 


14.78513590 


60.50336525 


-0.006867535404 


-0.003344481605 


27 


4.374333149 


15.96714209 


65.50388509 


-0.006104769071 


-0.002857142857 


29 


4.707718646 


17.14958715 


70.50433380 


-0.005481366771 


-0.002469135802 


31 


5.041096602 


18.33236558 


75.50472499 


-0.004963889096 


-0.002155172414 


33 


5.374468538 


19.51542724 


80.50506904 


-0.004528540629 


-0.001897533207 


35 


5.707835592 


20.69872494 


85.50537395 


-0.004157982647 


-0.001683501683 


37 


6.041198635 


21.88222143 


90.50564603 


-0.003839317503 


-0.001503759398 


39 


6.374558339 


23.06588687 


95.50589030 


-0.003562775224 


-0.001351351351 


41 


6.707915237 


24.24969671 


100.5061108 


-0.003320834926 


-0.001221001221 


43 


7.041269754 


25.43363246 


105.5063109 


-0.003107623257 


-0.001108647450 


45 


7.374622231 


26.61767652 


110.5064932 


-0.002918493910 


-0.00101112234 


47 


7.707972951 


27.80181575 


115.5066600 


-0.002749728194 


-0.000925925926 


49 


8.041322143 


28.98603869 


120.5068132 


-0.002598318198 


-0.00085106383 


51 


8.374670000 


30.17033570 


125.5069545 


-0.002461807379 


-0.00078492936 


53 


8.708016684 


31.35469854 


130.5070851 


-0.002338171717 


-0.00072621641 


55 


9.041362332 


32.53912015 


135.5072062 


-0.002225729996 


-0.00067385445 


57 


9.374707059 


33.72359445 


140.5073188 


-0.002123075281 


-0.00062695925 


59 


9.708050964 


34.90811616 


145.5074238 


-0.002029022010 


-0.00058479532 


61 


10.041394134 


36.09268071 


150.5075219 


-0.001942564749 


-0.00054674686 



In Table El we tabulate for odd d satisfying 5 < d < 61, the values, x max and x m i n , 
of x where the derivative Q' d {x) equals 0, the values of Qd(x) at these points, which 
will be absolute maximum and minimum values, and the limiting value Qn m of Q(x) 
as i -> ±oo. 

When d is even, the functions Qd fall into almost the same pattern, except that 
they have an additional wiggle between | and | + 1. For example, a schematic graph 
°f Q2&(z) is given in the bottom half of Figure 13.71 Note that the graph is drawn 
wildly out of scale. Similarly to Q25, it has an absolute maximum at x — 3.87 and an 
absolute minimum at x = 58.003. But instead of decreasing steadily between these, 
it has an additional single local minimum and local maximum which occur between 
x = 12 and 13. Maple calculations strongly suggest that for all even d, the graph 
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of Qd will have a form similar to that in the bottom half of Figure 13.71 and that 
the positions and values of the maxima and minima will have patterns extremely 
similar to those for odd d in Table [31 We will not pursue those here, as we prefer to 
concentrate on the simpler situation when d is odd. 

The reader will immediately be struck by the pattern in Table El which seems 
especially striking for x m i n (d). We have extended these calculations through d — 151, 
using 80 digits of accuracy in Maple. Then for k = 3, . . . , 10, we have found the real 
numbers cq, . . . , Ck which satisfy 

k 

(3.8) x min (rf) = |rf-2 + ^ ^ d = 51, 61,..., 51 + 10*. 

i=o ^ ' 

Each Cj seems to stabilize as k increases in f)3.8p . Moreover, using the formula (13.81) 
for x min (d) derived using just a few values of d gives agreement with computed val- 
ues of aJmin(rf) for all odd values of d to an increasing number of decimal places as 
k increases. In addition, we have, with k = 10, Cq equals, to 19 decimal places, 
.0104166666666666666 ^ 1/96. 

Conjecture 3.9. For odd d > 5, there are numbers x max and x m i n such that Qd{x) 
is decreasing for x max < x < x m i n , and increasing elsewhere. There are real numbers 
Ci for i > 1 such that 

oo 

Ci 

The initial digits of a, ... } c 6 are -.176504629629629, .16562740498, .20004439, .291872, 
.3215, and .28. 



Xmin(d) \d 2 + + ^ 



The 3-digit repetend in c\ leads one to guess that c\ = — 305/12 3 and that the Q 
are all rational numbers. 

We have performed similar analyses for x max and Qd(^min)- The initial terms are 
apparently x max = | - | - 64( ]_ 1} and Q d {x min ) = - • The series for ^^(^min) 
seems to converge more slowly than the others. 

We wish to emphasize that we cannot prove that Qd{x) for odd d has a unique 
maximum and minimum. This is all based on Maple calculations obtained by setting 
its derivative equal to 0, where Q d is a ratio of two polynomials of degree d — 1. 
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Table 4. Evidence for Conjecture 13.11 
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00401357 


26 





00403620 


27 





004045174 


28 





004042664 


29 





0040313148 


30 





004012644 
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003988285 



Conjecture 13.11 is equivalent to saying that Q d (x) < .005 for x > \{d + y/d + 2). 
This latter statement would follow from Conjecture 13 .91 expanded to include a formula 
for Qd(xmm) and to include even values of d, together with a proof that Qd{\(d + 
yjd + 2)) < .005. Some justification for this conjecture is given by Table HI which 
also shows why we use .005. 

There is one value of Qd(%), occurring just before the crucial range x > d/2, for 
which the value of Qd{%) is easily determined. This is given in the following result, 
whose easy proof we omit. 
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Proposition 3.10. If d is odd, then Qdi^f-) = 1- If d is even, then Qd(^) — 

-2 



(d+l)(d-2) ' 

4. The general case (si and s 3 not necessarily equal) 

In this section, we present our analysis of the general case, which is similar to, but 
not nearly so thoroughly developed as, the case s± = S3 considered in the preceding 
sections. 

For arbitrary sj, s 2 , and s 3 , now let 

i 

Maple suggests that for any values of the Sj this f Sl ,s 2 ,s 3 is a unimodal function of 
N. If so, we can find the value of N at which / achieves a maximum by an analysis 
extremely similar to that employed in the case s± — S3. 

The formula for / is symmetric in s\ and S3. We write s± = x and S3 = x + A, 
A > 0. Let d = si + s 3 - N, and 

Generalizing (12.41) . which is the case A = 0, we have 
(4 1) Pa ^ - 2x2 " 2(rf ~ A)x + ^ d ± 1 Z 2A) + A ' 1 



PA,d d + 1 PaAx) 

The easy generalization of (I1.7P is 

fx,s 2 ,x+A(2x + A-d) > f x ,s a ,x+A{2x + A - d - 1) 

iff 

( AO \ ^ 2z 2 +2(A+l)z+A 2 +A d+2 | RaA x ) 

^■ 2 > 52 - d+i 2-+ix^)> 

where RA,d{x) is the remainder in (14. ip . If we assume this remainder is negligible, 
then imposing equality in (14. 2[) and recalling si = x and s 3 = x + A yields 



d = J2s\ + 2 Sl + 2s| + 2s 3 + (s 2 + \f - s 2 - § 



and 



(4.3) N = si + s 2 + s 3 + | - J2s\ + 2sx + 2s 2 + 2s 3 + (s 2 + |) 2 , 
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as nice a generalization of II. 101 and (11.21) as one could possibly desire. This yields 



as the most likely number of elements in the union, assuming remainder terms are 
negligible. More analysis of the remainder terms is required. 

We have seen that when si = S3, the remainder terms can apparently only affect 
the value of AT by 1. In Table [5] we present data when S3 = Si + 8, indicating rather 
good agreement. Here "actual N" is where the maximum actually occurs. 

Department of Mathematics, Lehigh University, Bethlehem, PA 18015, USA 
E-mail address: dmdl@lehigh.edu 




MAXIMIZING A COMBINATORIAL EXPRESSION 



21 



Table 5. Comparison of actual N and formula N 
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